Attorney Docket: 46107-0037 
VEHICLE YAW STABILITY SYSTEM AND METHOD 
BACKGROUND OF THE INVENTION 
[0001] 1. Technical Field 

[0002] The present invention relates to a yaw stability system for a vehicle and, more 
particularly, to a yaw stability system having a control unit using a sliding mode control 
technique to determine a control yaw moment. 
[0003] 2. Description 

[0004] Yaw Stability Control (YSC) systems have been in use in the automotive industry 
for a number of years to increase the stability of the vehicle and to enhance vehicle 
performance. In general, YSC corrects the under- steering and over-steering of the vehicle in a 
handling maneuver (e.g. lane change, slalom, etc.), particularly on a low friction surface. It also 
helps the driver maintain yaw stability of the vehicle in a severe handling maneuver. 
[0005] The effectiveness of YSC systems varies widely depending on the system design. 
YSC systems commonly seek to minimize tracking error between a desired vehicle yaw rate 
and an actual vehicle yaw rate by selectively actuating braking mechanisms associated with the 
wheels of the vehicle. When a braking mechanism is actuated, the brake exerts a torque on the 
wheel which in turn induces a vehicle yaw moment. Most YSC systems are based on control 
techniques that rely on empirical data and are heavily dependent on testing. Systems based on 
M on-off ' control techniques commonly fail to consider the magnitude of the tracking error other 
than to determine the desired braking torque. For example, the systems do not "over-actuate" a 
braking device if the magnitude of the tracking error exceeds a predetermined threshold or 
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boundary layer. This deficiency in existing YSC systems often leads to undesirable braking 
device chatter as the devices are repeatedly actuated. 

[0006] Moreover, conventional YSC systems commonly use hydraulically actuated 
friction based braking devices to induce control yaw moments. While these systems are 
generally suitable, they suffer from undesirably long response times and lack of smoothness 
during operation. The harsh operation of the braking systems induce undesirable noise, 
vibration, and harshness (NVH) during operation. 

[0007] As is discussed in detail in this application, one feature of the present invention is 
the use of electromagnetic retarders, preferably eddy current machines, as YSC braking 
devices. While electromagnetic retarders have been used in braking systems for commercial 
trucks for many years, these retarders are generally not used in YSC systems for a number of 
reasons, including difficulty in accurately modeling the torque characteristics of the retarder. 
One modeling consideration of particular interest in YSC systems is the ability to obtain an 
accurate estimation of the retarding torque generated by an electromagnetic retarder. Accurate 
torque estimation is important for providing consistent performance. One conventional 
estimation technique requires an initial estimation of armature temperature which is then used 
in the torque calculation. Others have estimated electromagnetic retarder braking torques using 
predetermined look up tables of torque versus peak voltage between the retarder poles at 
various rotor speeds. Yet others have modeled eddy current brakes as a function of excitation 
current and rotor speed. However, each of the aforementioned techniques suffer from 
inaccuracies, assumptions that are not appropriate for many operating conditions, and/or 
computational intensity. 
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[0008] Thus, a need exists for an accurate and systematic YSC approach that minimizes or 
eliminates undesirable chattering, reduces dependence on empirical data and testing, improves 
response time, and minimizes NVH. 

SUMMARY OF THE INVENTION 
[0009] The present invention is generally directed to a robust YSC method that improves 
vehicle yaw stability performance. Another aspect of this invention is a vehicle yaw stability 
system that utilizes electromagnetic devices to generate a control yaw moment. The YSC 
system of the present invention is based on a sliding mode control theory and, in the illustrated 
embodiment, utilizes eddy current electromagnetic retarders. The rapid response from the eddy 
current retarders is aided by the sliding mode control law and the smooth control of the 
electromagnetic machines provides improved NVH relative to hydraulically actuated friction 
braking devices. Further, the robustness in the controller maintains good performance even in 
the presence of external disturbances. 

[0010] In one embodiment, the invention is directed to a yaw stability system for a vehicle 
having a plurality of wheels each with a torque element. The yaw stability system includes a 
yaw rate sensor, a plurality of braking devices, and a control unit. The control unit 
communicates with the yaw rate sensor and is configured to identify a desired yaw rate, 
determine a yaw rate tracking error, determine a yaw condition of the vehicle based on the 
vehicle yaw rate, determine a control yaw moment to minimize the yaw rate tracking error, 
select one or more of the plurality of braking devices based on the yaw condition, and 
communicate a control command to one or more of the selected braking devices to induce the 
control yaw moment. 
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[0011] In another embodiment, the present invention is directed to a method for 
controlling yaw in a vehicle. The method includes determining a vehicle yaw rate and desired 
yaw rate, calculating a yaw rate error, determining a control yaw moment using a sliding mode 
control law based on a lumped mass vehicle model, selecting one of the braking devices based 
on a vehicle yaw condition, determining a control command based on the control yaw moment, 
and communicating the control command to one of the selected braking devices. 
[0012] In still another embodiment, the present invention is directed to a method of 
estimating the retarding torque of an electromagnetic retarder. The method includes modeling 
the estimated retarding torque based on a quadratic function of a control element velocity and a 
magnitude of a retarder excitation current. The method further includes estimating coefficient 
functions of the quadratic relationship from steady state test data at various rotor speeds. 
[0013] In a further embodiment, the invention is directed to an electromagnetic retarder 
assembly including an electromagnetic retarder, a sensor, and a controller. The 
electromagnetic retarder has a stator with conductive windings and a rotor rotatable relative to 
the stator. The sensor senses the rotational speed of the rotor and generates signals indicative 
of the rotational speed. The controller communicates with the sensor to receive the signals and 
is configured to communicate an excitation current to the electromagnetic retarder. The 
controller is also configured to estimate the retarding torque of the retarder based on a quadratic 
function of the velocity of the rotor and the magnitude of the excitation current. 
[0014] Further scope of applicability of the present invention will become apparent from 
the following detailed description, claims, and drawings. However, it should be understood 
that the detailed description and specific examples, while indicating preferred embodiments of 
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the invention, are given by way of illustration only, since various changes and modifications 
within the spirit and scope of the invention will become apparent to those skilled in the art. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0015] The present invention will become more fully understood from the detailed 
description given here below, the appended claims, and the accompanying drawings in which: 
[0016] Figure 1(a) is a schematic illustration of a vehicle having a YSC system according 
to the present invention and representing vehicle yaw dynamics in a non-braking condition; 
[0017] Figure 1(b) is a schematic illustration similar to Figure 1(a) illustrating control yaw 
moment generation via braking forces; 

[0018] Figure 2 is a schematic overview of the control unit control strategy; 

[0019] Figure 3 is a flow chart illustrating a suitable control strategy for a control unit of 

the present invention; 

[0020] Figure 4 is a sectional view of a representative eddy current machine for the torque 
actuating devices shown in Figure 1 ; 

[0021] Figure 5 is a graph showing the retarding torque versus current characteristics of 
the eddy current machine at various rotor speeds between 100 and 1000 RPM; 
[0022] Figure 6 is a graph showing the torque versus rotor speed characteristics of the 
eddy current machine at constant excitation current; 

[0023] Figure 7 includes figures 7(a), 7(b), and 7(c) showing the plots of coefficient 
functions fo,fi, and ^ for a variety of rotor speeds as identified from the first set of least square 
fit for estimating the parameters aoo, a>ou • • •> m & a 22\ 

[0024] Figure 8 is an overlay plot of estimated torque and measured torque for the front 
wheel retarders; and 
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[0025] Figure 9 is an overlay plot of estimated torque and measured torque for the rear 
wheel retarders. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0026] Figure 1(a) schematically illustrates a vehicle 10 having a plurality of wheels 12 
rotating with shafts 14. The vehicle includes a YSC system 16 illustrated to include braking 
devices 18a-18d associated with each of the wheels, a yaw rate sensor 20, and a control unit 22. 
During operation, any change in the direction of the vehicle 10 generates a yaw rate about the 
vehicle's center of gravity (CG) 24. While a certain yaw rate is desired for proper vehicle 
turning, the vehicle operator may place the vehicle in oversteer or understeer conditions where 
the vehicle yaw rate differs from the desired yaw rate. The control unit 22 seeks to identify and 
correct the tracking error between the measured vehicle yaw rate and the desired yaw rate. 
[0027] The YSC system 16 uses a systematic approach to control vehicle yaw rates. This 
systematic approach, based in part on using a sliding mode control law in combination with a 
lumped mass vehicle model to determine a control yaw moment, improves the system response 
time, enhances vehicle yaw rate tracking to the desired yaw rate, and improves controller 
robustness and stability when compared to conventional systems. The YSC system described 
in detail herein uses electromagnetic retarders, particularly though not necessarily eddy current 
machines, as the braking devices 18. While the use of electromagnetic retarders provides 
numerous operational benefits over conventional systems it should be recognized that the 
control strategy may be used with YSC systems incorporating other braking mechanisms, 
including friction brakes. 

[0028] An electromagnetic retarder or brake 30 suitable for use as a braking device is 
illustrated in Figure 4. Retarders of the type illustrated in Figure 4 and described herein are 
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generally known in the art and follow the basic principles of electromagnetic induction. In 
general, the retarder 30 (illustrated as an eddy current machine) has an iron core 32, a stator 34, 
conductive windings 35, and a rotor 36 fixed to rotate with the wheel shaft (not shown). 
Providing the stator windings 35 with an excitation current induces an eddy current in the stator 
and a retarding force which acts on the rotor 36. A conventional sensor monitors the rotational 
speed of the rotor 36 and communicates a rotational velocity signal to the control unit 22. The 
electromagnetic retarders may be used as secondary retarders for each vehicle wheel or as the 
primary braking system for the vehicle. 

[0029] The operation of the control unit 22 will now be described with reference to the 
flow charts illustrated in Figures 2 and 3. As an overview, the control unit 22 receives input 
from conventionally configured sensors or stored data and determines and communicates a 
current command to one or more of the electromagnetic retarders 30 to minimize yaw rate 
tracking error. The control unit 22 is generally configured in a conventional manner to carry 
out the tasks described herein. Thus, it will be appreciated that the control unit may take a 
variety of forms without departing from the scope of the present invention. By way of 
example, it is contemplated that the control unit includes a microprocessor with a calculation 
module, a memory or data structure for storing and/or retrieving data, as well as appropriate 
input and output circuits for receiving the various input signals and communicating control 
commands to YSC system components such as the retarders 30. 

[0030] For ease of explanation, the control strategy is separated into four stages: (1) 
determining the yaw condition (i.e., understeer or oversteer) of the vehicle (Step 40 in Figure 2) 
such as by comparing the vehicle yaw rate determined from the yaw rate sensor 20 to a desired 
yaw rate; (2) determining a control yaw moment (M z ) if the yaw rate error (i.e., the difference 
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between the desired and measured yaw rates) exceeds a threshold (Step 42); (3) determining a 
torque command for the braking devices 18 (Step 44) to correct the error; and (4) calculating an 
appropriate control command for one or more of the braking devices (Step 46). While each of 
these steps are described below and illustrated in Figures 2 and 3 in the sequence set forth 
above and using specific relationships and equations, it should be appreciated that the sequence 
of the steps and the manner in which the determinations are made may be varied without 
departing from the scope of the invention defined by the appended claims. By way of example 
rather than limitation, it is noted that the determination of the control yaw moment (Step 42) is 
described as being performed using the relationship set forth below in Equation (8) below and 
after the determination of the yaw condition (Step 40). However, as the described manner of 
determining the control yaw moment is independent of the yaw condition, the control yaw 
moment may be determined before the determination of the yaw condition. Moreover, while 
the relationships used to determine the control yaw moment provide operational benefits 
related, for example, to calculation efficiencies and accuracy, other methods for calculating the 
control yaw moment (M z ) may be used. 

[0031] The YSC system has access to various conventional sensors or stored data to 
facilitate the performance of the calculations discussed herein. For example, the control unit 22 
is shown in Figures 1 and 3 to communicate with a steering wheel angle sensor (50), the 
vehicle yaw rate sensor (20), and wheel speed sensors 52 as well as to have access to estimated 
road friction coefficients, wheel patch normal forces, and slip angles. A variety of suitable 
techniques for measuring or estimating these parameters are generally known in the art and 
widely published in available literature. One skilled in the art will also recognize that 
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conventional mechanisms (such as Kalman filter, fuzzy logic, or observer) may be used to 
tailor the parameter estimation for a particular application. 
[0032] L Determination of yaw condition 

[0033] As noted above, the control unit 22 of the YSC system seeks to minimize the 
tracking error between the desired vehicle yaw rate and the measured yaw rate. The tracking 
error is represented below in equation (1) as the sliding or yaw rate error surface (S): 

S = r des (t)-r(t) (1) 

where 

r des (0 = Desired vehicle yaw rate at time t; and 

r(t) = Vehicle yaw rate at time t. 
[0034] In the preferred embodiment of the present invention, the desired vehicle yaw rate 
is determined for a number of vehicle operating conditions and stored in a data structure (such 
as via a look-up table) with which the control unit 22 communicates. The desired values for the 
yaw rate are preferably based on experimental measurements performed prior to 
implementation of the system in the vehicle. Those skilled in the art will appreciate that a 
variety of techniques may be used to determine and store desired vehicle yaw rate values for 
access by the controller. The vehicle yaw rate at any given time is preferably obtained from the 
yaw rate sensor 20 shown in Figure 1. 

[0035] In Step 60 (Figures 2 and 3), the control unit 22 determines whether the tracking 
error exceeds a threshold. This threshold may be determined in any number of ways known in 
the art. In general, the threshold is determined so as to minimize calculation complexity 
without negatively impacting vehicle performance. Factors impacting the selection of an 
appropriate threshold include the occupants' perception of understeer/oversteer conditions, 
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vehicle stability limits, the value of the measured yaw rate, road friction coefficients, etc. 
Accordingly, it is anticipated that the control unit 22 will communicate with a data structure 
having a look-up table identifying the threshold for a number of vehicle operating conditions. 
[0036] As is indicated in Step 60 of Figure 3, if the measured vehicle yaw rate does not 
exceed the threshold, the control unit 22 proceeds no further through the control strategy and 
performs the threshold comparison at predetermined intervals. In the event the yaw rate 
exceeds the threshold, the control unit 22 proceeds to Step 40 to determine the yaw condition of 
the vehicle. With input from the vehicle yaw rate sensor 20, steering wheel angle sensor 50, 
and information regarding the desired yaw rate, the control unit 22 determines the direction of 
the steering wheel angle (Step 62) and the vehicle yaw condition (Steps 64 and 66). There are 
two situations that accompany yaw instability - understeering and oversteering. In an 
understeer condition the absolute value of the vehicle yaw rate, r, is smaller than the absolute 
value of desired vehicle yaw rate, rd es . Conversely, in an oversteer condition, the absolute value 
of the vehicle yaw rate, r, is larger than the absolute value of desired vehicle yaw rate, rd es . 
Based on this relationship, the control unit 22 determines in Steps 64 and 66 whether the yaw 
instability is due to an understeer or oversteer condition. 
[0037] 2. Determination of control yaw moment. 

[0038] The determination of the control yaw moment, M z , is represented by Step 42 in 
Figures 2 and 3. For clarity, the determination of the control yaw moment in Step 42 is shown 
as separate calculation Steps 70, 72, 74, and 76 for each of the yaw conditions determined in 
Step 40. However, as indicated in Figure 3 and described herein, the control yaw moment in 
each instance is preferably calculated based upon Equation (8) below. 
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[0039] The control unit 22 may use a variety of vehicle yaw dynamic relationships 
generally known in the art to determine the control yaw moment. In the embodiment 
specifically described herein, the control unit 22 uses a sliding mode technique based on 
lumped mass vehicle model. Unlike conventional YSC systems (such as those based on 
Proportional Integral and Derivative Control (PIDC)), use of a sliding mode technique based on 
a lumped mass vehicle model for the YSC system provides fast response and robustness with 
respect to external disturbances. As fundamentals of the sliding mode control law are generally 
known in the art, all details of the technique are not described in detail herein. However, 
Equations (2) through (8) and the related descriptions illustrate the technique in a manner that 
is sufficient for use with the invention. Using the lumped mass model, the vehicle yaw 
dynamics can be obtained as follows: 

dr 

*zz— = *(F' yFL cos 5 \ + F \FL sin S \ + F 'yFR cos 8 2 + F xFR sin 5 2 > " b ( F ' yRR +F * yRL > (2) 
+ c(F' yFL sin S x - F xFL cos - F< xRL ) - d(F< yfR sin S 2 - F xFR cos S 2 - F< xRR ) + M z 

where 

Izz = Vehicle yaw inertia; 
M z = Control yaw moment; 

F*xfl, FVfl, Fxfr, F'yFR, F x rl, F'yRL, F' y RR, F' y RR, = Tire contact patch forces in x- and y- 

directions in a non-braking condition as illustrated in Figure 1(a); 
8i, 62 = Road wheel angle for the front wheels; and 

a, b,c,d = Contact patch locations from the vehicle center of gravity ("CG") 
[0040] In equation (2), the road wheel angle for the front wheels (12a, 12b) may be 
determined through the use of conventional steering wheel angle sensor 50 or road wheel angle 
sensor 52. The contact patch locations relative to the vehicle CG are illustrated in Figure 1(a), 
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are known, and may be stored for use by the control unit in any conventional manner such as in 
a look-up table. To simplify the calculation of equation (2), the control strategy assumes that 
the road wheel angle for the front left tire 12b is equal to the road wheel angle for the front 
right tire 12a and that the wheel separation is the same for the front and rear wheels. 
[0041] Rewriting equation (2), we obtain, 

IJ = 4TyFL +F' yFR )cos6 + (F< xFL + F' xFR )smS\-b(F' yRR +F' y J+ 

c{f\ fl sin 8 - F\ FL cos 8 - F' xRL ) - d(F\ FR sin 8 - F\ FR cos 8 - F' xRR )+ M 2 

[0042] As is shown in equation (4) below, control unit 22 operates on the following further 
assumptions: (1) that the normal force (F zX ) on the left and right side of the vehicle is same 
{i.e., the normal force on the front left contact patch is the same as that on front right contact 
patch, etc.); (2) that the rolling resistance coefficients (r/ij) dictate longitudinal forces at tire 
contact patches; and (3) that the lateral friction forces vary linearly with the slip angle (a) for 
each wheel. 

F zfl~ F z fr = F zF \F zRl = F zRR = F zR 

F xFL = iFL^zF > F xFR = 7 7 FR^zF ^ xRl ~ Vrl^zR > ^ xRR = ^RR^zR ^ 
a FL =CC FR =a F > a RL ~ a RR =CX R 

F yFL ~ C F L ctf ; F yFR = C FR a F \F* yRL = C RL a R \F yRR = C RR a R 

where 

Cfl, Cfr, Crl, & Crr are the cornering coefficients from a two track vehicle model 
which may be determined based on experimental data in a manner generally 
known in the art. 

ocfu ocfr, a RL> & ocrr are slip angles associated with each wheel determined in a 
conventional manner. 
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Vfl, Vfr, ??rl 9 & ??rr are the rolling resistance coefficients associated with each road-tire 
contact patch determined in a conventional manner. 
The simplification set forth in Equation (4) provides the following yaw dynamics equation for 
determining the rate of change of the vehicle yaw rate ( r ): 

r = j- a[{C FL + C FR )a p cosS + (f, pL + ^ sin J] - b(C RR + C^)<^ + (c * - * C^)^ sin* - 

zz P) 

cosS + lRL F zR > + d ^F^ zF cosS + iRlFJ + M * 
[0043] As noted above, the control unit 22 uses a sliding mode control law which, in 
general, generates the following relationship between the sliding surface and its gradient (S ): 

S = -7jSAT - 
\<PJ 

which, when substituting the relationship of Equation (1), provides: 



r des -r = -tjSAT(^— (6) 



V 



where 

S = is the sliding or yaw rate error surface as set forth in Equation (1); 
r| = Convergence factor; 

<|> = Boundary layer thickness of control command determined; and 

SAT = Saturation function 
The convergence factor and boundary layer thickness in Equation (6) each impact the 
performance of the system, particularly with respect to response time and chatter. As control 
within the boundary layer is smooth and stable, it is generally desirable to select parameters 
that rapidly bring the control within the boundary layer and maintain the control within this 
layer during operation. Suitable values for these parameters may be determined through 
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tuning, such as via computer assimilation modeling and/or system prototype testing. The 
saturation function is a mathematical function that sets the controller output at +/- 1 if the 
output is outside a range where the output behaves linearly. By way of example, it has been 
determined that a convergence factor within the range of about 50 to about 100 and a boundary 
layer thickness between about 5 and 15 degrees per second are suitable for a yaw stability 
control of a light duty truck. Notwithstanding this illustrative example, those skilled in the art 
will appreciate that the parameters may vary based on the application and the desired 
performance. Substituting the relationships of Equation (5) into Equation (6) yields: 

M z = tJdes ~ i a [( C FL + C FR ) a F cosS + ^FL + ^FR^ F zF sin ^ " b ^ C RR + C R^ a R + 

r d -r (7) 

(c*C FL -d* C FR )a p sin S - c^F^ cos* + rj^ ) + d^F^ cos* + ^FJ] + IjfiA Vr*^) 

If the rolling resistance forces at the tire contact patches are neglected, the following control 
yaw moment (Mz) is obtained from Equation (7): 

K =I J<ies -HC FL +C FR )a p cosS-b(C RR +C RL )a R +(c*C FL -d*C FR }a p vnS^I^TjSAT^^) (8) 

[0044] Each of the parameters set forth in Equations (7) and (8) are known or determinable 
by the control unit 22. In the illustrated example: 

Tire road rolling resistance coefficients (ijij) for each contact patch are estimated by 

conventional techniques; 
Normal force (F zF ) are assumed to be the same on the left and right side of the vehicle 

and determined by conventional means; 
Cornering coefficients {C FL > Cfr, Crl, and Crr) are preferably stored in the controller 
database after determination by conventional means such as a two track vehicle 
model; 
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Slip angles associated with each wheel (otp and aft) are assumed to be the same for each 

wheel on the front and rear axles and are determined by conventional means; 
The vehicle yaw rate (r) is preferably provided by the yaw rate sensor 20; 
The desired vehicle yaw rate (r des ) is determined by experimental data or data from 

previous developments in a manner generally known in the art; 
The rate of change of the vehicle yaw rate (r des ) is calculated from two successive 

measurements of the desired vehicle yaw rate (r^); and 
The saturation function (SAT) represents a value of +/- 1 when the absolute value of the 

argument of the function is greater than unity, otherwise it assumes the value of 

the argument itself 

Notwithstanding this representative illustration of the preferred manner of obtaining 
measurements or values for calculating the control yaw moment from Equation (8), those 
skilled in the art will appreciate that a variety of measurement or estimation techniques may be 
used with the present invention. 

[0045] It is further noted that the term I zz represents the vehicle mass moment of inertia 
about a vertical axis through the vehicle center of gravity 24. While this term is generally 
constant for the vehicle, its value may vary slightly with varying vehicle load. For simplicity, 
the control unit 22 of the present invention assumes that the term is constant and stores the term 
in its data structure. The term F zF represents the normal force on the front contact patch and is 
determined by dynamic load distribution. Thus, with the other parameters of Equations (7) or 
(8) being known or determinable in the manner discussed above, the control unit 22 can 
determine the control yaw moment (M z ). That is, the yaw moment to induce in order to 
minimize the tracking error relative to the desired yaw rate. 
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[0046] It is noted that the above method of determining the control yaw moment provides 
advantages not previously recognized in the art. For example, the derivative of the desired yaw 
rate adds a dynamic term to the control command in order to improve responsiveness. Further, 
to simplify the control strategy and reduce response time, Equations (7) and (8) includes the 
saturation function. Notwithstanding the benefits of the above described strategy for 
determining the control yaw moment, it will be appreciated that other strategies and methods 
maybe used without departing from the scope of the invention defined by the appended claims. 
[0047] 3. Determination of Torque Command 

[0048] With the vehicle yaw condition and control yaw moment determined as described 
above, we now turn to the determination of the torque command communicated to the braking 
devices in order to induce the desired control yaw moment. Equations (10) through (19) 
represent the control law for the yaw control system proposed of the present invention. Step 44 
represents the determination of the torque or control command based on the control yaw 
moment, M z , calculated above. That is, the control unit 22 determines which of the braking 
devices 18a-18d to engage and the magnitude of the braking forces required to generate the 
control yaw moment. The control unit then communicates an appropriate torque or control 
command to one or more of the braking devices. 

[0049] As noted above, the YSC system may include different types of braking devices for 
generating the control yaw moment. The following description is provided with reference to 
the illustrated embodiment of the YSC system using electromagnetic retarders 30, preferably 
eddy current machines, operatively associated with each of the vehicle wheels 12 and, more 
particularly, with the rotors. The control yaw moment is generated by selectively energizing 
one or more of the electromagnetic retarders 30. 
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[0050] In an understeer condition, the control yaw moment is generated by applying a 

braking force to one or more of the inner wheels (e.g., left side wheels 12b and/or 12d during a 

left hand turn) whereas in an oversteer condition the control yaw moment is generated by 

applying a braking force to one or more of the outer wheels (e.g., the right side wheels 12a 

and/or 12c during a left hand turn). In any of these two vehicle dynamic conditions, either both 

wheels or one wheel (on one side) can be braked to generate the desired control yaw moment. 

However, to minimize control complexity it is desirable to brake only one wheel to generate the 

control yaw moment. In the case of braking only one wheel, it is more effective to brake the 

front wheel in an oversteer condition and the rear wheel in an understeer condition. 

[0051] With the above in mind, and assuming counterclockwise rotation positive and that 

the slip angles are small compared to road wheel angles, the relationship between the control 

yaw moment and brake force is, with reference to Figure 1(b), as follows: 

M \ = cF xFL cos S - aF xFL sin S - dF xFR cos 5 - aF xFR sin 8 + cF^ - dF RR 

' (9) 

where F x fl, F x fr, Frl, Frr are the braking friction forces at the tire contact patches. 

[0052] Equation (9) may be solved for instances where a braking force is applied to any of 

the four wheels. As is generally known, the torque induced by any braking device in the YSC 

system is a function of the rotational speed of the control element (e.g., rotor). Therefore, 

while it is generally desirable to brake only a single wheel, a single braking device may not be 

capable of generating the desired control yaw moment. As to the illustrated embodiment using 

electromagnetic brakes, the relationship between the induced torque and control element speed 

may, in certain situations, cause a retarder to saturate. Saturation of electromagnetic retarders 

occurs when a maximum control current is provided to the retarder. In these instances, both 
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front and rear wheel braking devices on a given vehicle side can be used to generate the desired 
yaw moment. 

[0053] A. Single Wheel Braking 

[0054] Equations (10) - (13) indicate the torque generated by braking one of the wheels 
12a-12d in response to an understeer or oversteer condition based on the relationship set forth 
in Equation (9). The torques set forth in each of Equations (10)-(13) are a function of the 
wheel rolling radius, represented by "R". 

[0055] In an understeer condition, the torque generated by braking the left rear wheel 12d 
during a left hand turn (i.e., vehicle turning counterclockwise) may be determined in Step 80 
(FIG. 3) by Equation (10): 

M 2 =cF RL = c ^. 

R 

T tRL = —M 2 

c 

(10) 

[0056] Similarly, in an understeer condition, the torque generated by braking the right rear 
wheel 12c during a right hand turn (i.e., vehicle turning clockwise) may be determined in Step 
82 by Equation (11): 



M z =dF RR =d- 



T bRR =^ M * 



bRR 

R 



(11) 
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[0057] In an oversteer condition, the torque generated by braking the right front wheel 12b 
during a left hand turn (i.e., vehicle turning counterclockwise) may be determined in Step 84 by 
Equation (12): 



M z ~ (d cos 8 -a sin S)F xFR = (d cos 8 - a sin 8) 



bFR 

R 



T bFR = M 2 (12) 

bFR (d cos 8- a sin 8) 2 

[0058] Finally, in an oversteer condition, the torque generated by braking the left front 
wheel 12b during a right hand turn (i.e., vehicle turning clockwise) may be determined in Step 
86 by Equation (13): 

T 

M z = (ccosS -a sin S)F xFL = (ccos£-a sin£)-^- 

R 

(c cos 8 - asm 5) 

[0059] Each of the above equations (10) - (13) indicate the torque required from the 
respective braking devices 18a-18d to induce the control yaw moment (M z ) the braking of a 
single wheel. For completeness, it is noted that if the denominator of Equations (12) or (13) are 
zero, the denominator is set to a small value, indicated by TOL, to prevent overflow conditions 
in the control unit. 

[0060] B. Torque Saturation 

[0061] Steps 90, 92, 94, and 96 in Figure 3 represent the comparison of the required torque 
from a single braking device (e.g., TbRL in Step 90) to the estimated torque that the braking 
device is capable of generating (e.g., T estRL ). The torque generated by the electromagnetic 
retarders 30 is a function of rotor speed and the current supplied to each device. This 
relationship is illustrated by the curves shown in Figures 5 and 6. More particularly, Figure 5 
shows the retarding torque versus current characteristics of the eddy current machine at various 
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rotor speeds between 100 and 1000 RPM. Figure 6 shows the torque versus rotor speed 
characteristics of these machines at constant excitation current. As indicated in Figure 3, if the 
required torque is less than the estimated torque, the control unit 22 proceeds to Step 46 and 
calculates the current command according to Equation (19) below. However, if the required 
torque is greater than the estimated capacity of the actuator, actuator saturation prevents the 
inducement of the required control yaw moment. 

[0062] In identifying retarder saturation, the control unit 22 uses a torque estimation 
algorithm to determine the estimated torque (T est ) which an eddy current machine can generate 
based upon the rotor speed (o>) and excitation current (i). Previous modeling of 
electromagnetic retarders, particularly eddy current retarders, have assumed that the braking 
force or torque varies linearly with speed. In the illustrated embodiment of the present control 
strategy, the braking torque is more accurately modeled through a quadratic function of rotor 
speed and excitation current. The control unit 22 is preferably an open loop controller 
providing a current optimal torque without the need for a current feedback signal from the 
retarder thereby eliminating the need for a current sensor. Further, the invention contemplates 
use of a computationally efficient parametric model control unit to make the YSC system 
suitable for implementation as production hardware. For example, a parametric model based 
controller requires less memory in the microprocessor as compared with a look-up table based 
controller. Hence both the computational and hardware cost can be reduced with a parametric 
model based controller. 

[0063] The steady state torque model for the eddy current machines is represented by the 
following Equation (14). 

7 - =/.(*) + / I (»)*i + / 2 (ffl)*/ a (14) 
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where 

T = retarding torque 

i = retarder feedback current 

f(co) = a i0 +a n o + a i2 a) 2 

with 

ax) , an , aa = identified parameters 
co = rotor speed 

The parameters ay (where i = 0, 1, 2 and j = 0,1,2) in Equation (14) are functions of the rotor 
speed and are estimated from steady state test data. The solid lines 110, 112 and 1 14 in Figures 
7(a), 7(b), and 7(c), respectively, show plots of coefficient functions fo,fj, and f 2 for each rotor 
speed as identified through a least square estimate based on the steady state test data. The 
solid line plots 110, 112, and 114 of f 0 , //, and f 2 are used for the initial estimate of the 
parameters a 0 o, a 0 j, and a 22 . The coefficient functions fo, fi, and f 2 are then re-calculated 
and plotted for each rotor speed based on the initial estimates of parameters ay. The plots of 
these functions are shown by dotted lines 116, 118, and 120 in Figures 7(a)-(c). The 
parameters ay for each of the coefficient functions fo, fu and f 2 are then re-estimated through 
another least square type algorithm. A variety of conventional least square fitting techniques, 
including the MATLAB function PLYFIT, may be used for each of the above discussed least 
square estimations. It is apparent that the second stage fit of the coefficient function /o,//, and f 2 
(shown by dotted lines 116, 118, and 120) match reasonably well with the first stage estimate 
of the coefficient functions directly from experimental data. Accordingly, Figures 7(a)-(c) 
graphically illustrate that the torque produced by the retarder 30 is accurately modeled using a 
quadratic function of the excitation current and the rotor speed. 
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[0064] If an eddy current machine is used as the primary vehicle braking system, the 
desired torque in a brake-by-wire system may be tracked in order to meet the driver's brake feel 
requirements (pedal displacement/force versus vehicle deceleration). It is also desirable, 
whether the eddy current machine is used as a primary vehicle braking system or as a 
secondary system in combination with conventional hydraulically actuated friction brakes, to 
provide accurate torque control in an anti-lock brake system, traction control system, or a 
vehicle stability control system. The estimated torque produced by the electromagnetic retarder 
may be calculated based on the retarder torque model discussed in detail below, particularly 
those sections referring to Equation (14). The accuracy of this estimation is evidenced by the 
overlayed plots of estimated torque and measured torque for the front and rear wheel retarders 
(Figures 8 and 9, respectively). For completeness it is noted that the measured wheel torque is 
the net torque at the wheel indicated by a torque sensor. The illustrative tests were designed to 
minimize the effects of driveline torque and inertia on the measured wheel torque. The test 
conditions included performing all braking tests in neutral gear, limiting tests to straight line 
braking with only eddy current machines as the retarding device, and dry asphalt road 
conditions. 

[0065] Using Equation (14), the control unit estimates the torque that an electromagnetic 
device is capable of generating at a given rotor speed and current. As noted above, if the 
estimated torque for the desired wheel is less than the torque required pursuant to (10) through 
(13), then the control unit 22 determines the amount of torque to be generated by the other 
wheel on the same side of the vehicle as described below. 
[0066] C. Multiple Wheels 
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[0067] When the braking device associated with the desired single wheel is incapable of 
generating the desired torque prior to saturation, both front and rear wheel braking devices can 
be used to generate the requested torque. The control unit 22 is configured to identify when 
retarder saturation occurs (Steps 90-96) and to generate current commands to the retarders 
associated with each wheel on a side of the vehicle when the control yaw moment would 
saturate a single actuator. Equations (15) - (18) below may be used to determine the amount of 
torque to be generated by each wheel. First, the torque to be generated by the wheel that is 
most effective at correcting the tracking error (e.g., the rear left wheel for counterclockwise 
understeering) is set at the estimated maximum or saturation level (Test) calculated above. As 
illustrated in Steps 100-106 and described below by Equations (15) - (18) the control unit then 
determines the amount of torque to be generated by the next most effective wheel (e.g., the left 
front wheel in a counterclockwise understeering condition). 

[0068] With the above in mind, Equation (15) sets forth an equation for the braking of the 
front and rear left wheels for a counterclockwise understeered condition. 



T T 



M 2 = (ccos8-asmS)F xFL +cF RL = (ccos£-asin S)^^ + c-^- 

R R 

JfTbRL > T es tK L ^ nT bRL = T cs,RL 

_ RM z — cT eslR L 

bFL ' 



(15) 



(ccos 8 -asm 8) 

Similarly, Equation (16) sets forth an equation for determining the front and rear right wheel 
braking torques for an understeered vehicle turning clockwise. 

M z = -{d cos 8 + asm S)F xFR — dF RR = -(</cos 8 ^ asmS)^-d^ bRR 



R R 

V^bRR > ^ estRR > thCH T = T £S(RR 

RM +dT 

7 - l estRR 

bFR (dcos8 + asmS) 
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Further, Equation (17) may be used for determining the brake torque for the front and rear right 
wheels when the vehicle is oversteered and turned counterclockwise. 



M z = -(dcosS + asinS)F xFR -dF RR =-(rfcos£ + asin S)^- d^- 

R R 

IJTbFR > T eslFR ,thenT hFR =T eslFR 

RM z +(dcosS + a sin S)T eslFR 
d 



(17) 



Finally, Equation (18) may be used to determine the braking torque for the front and rear left 
wheels in a clockwise oversteered condition. 

M z = (c cos 5 - a sin 8 )F xFL + cF RL = (c cos 8 - a sin 8) + c ? hRL 



R R 

WbFL > T es<FL> thenT bFL ~~ ^ estFL 

RM z - (c cos 8 - a sin S)T eslFL 
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Again, for completeness, it is noted that if the denominator of Equations (15) or (16) in Steps 

100 or 102 is equal to zero, the denominator is set to TOL. 
[0069] 4. Generating current command 

[0070] Once the torque command has been calculated for each of the selected 
electromagnetic retarders based on the above equations, the control unit 22 communicates a 
current command to the respective retarder, e.g., eddy current machine, in (Step 46) for a given 
wheel speed (assuming that the retarder is not saturated) and is represented by Equation (19): 



/?(*>) 

where 

Ixy = Current command to FL, FR, RL, or RR eddy current machine 
Tbxv = Desired torque for the FL, FR, RL, or RR eddy current machine 
fo y j f\ Xy = Speed dependent retarder parameters 
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[0071] Equation (19) is derived from the relationship set forth in Equation (14), which 
established the quadratic equation as a function of current and rotor speed, with the assumption 
that the measured current is the same as the commanded current, which is particularly 
applicable for relatively short time constants as in the illustrated embodiment. More 
particularly, given the desired torque and measured rotor speed, an analytical quadratic 
equation solver may be used to obtain the current command from Equation (14) as follows: 



_ -fM+Jf? (a>)-4/ 2 (Q>)(Mco)-tZ) 



r _ /.(«) 

1 cmd ~ 



2/ 2 (») 

( 



(20) 



V 



l 4/ 2 M(/ 0 (^)-r^) 

A 2 (CO) 



2f 2 (o>) 

If the second term in the square root symbol is much smaller than unity, the term under square 
can be expanded in a Taylor series form as follows: 

(l + xh =1 + + 

v ' 2 2! 4 

Using the above expansion, Equation (20) is modified and re-written as follows: 

_ Mco) ( 2/ 2 ( fl Q(/ 0 (< P )-r fl .) ) 

2f 2 (co){ f x (a)) J 

If the second term in the square bracket of Equation (20) is much smaller than unity, the higher 
order terms can be omitted. Simplifying the above equation yields Equation (19) which is the 
proposed open loop control law for an eddy current brake system. The coefficient function fo 
and fj are estimated through a parameter estimation scheme as mentioned above. Appropriate 
checks should be performed for coefficient function fj in order to prevent any possible overflow 
condition. 
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f x {a>) 

[0072] From the above description, it will be appreciated that the YSC system and method 
of the present invention provides numerous advantages over conventional systems and 
techniques. By way of example, the present invention utilizes an improved control yaw 
moment calculation methodology, based on a sliding mode control law based and a lumped 
mass vehicle model, to minimize tracking error between a desired vehicle yaw rate and an 
actual vehicle yaw rate. This method provides smoother, more stable, and more robust control 
so as to minimize undesirable braking device chatter. While the YSC system of the present 
invention may be used with a variety of braking devices, the disclosed embodiment particularly 
relates to the use of electromagnetic braking devices, and more particularly to eddy current 
machines, so as to provide improved response times and undesirable noise, vibration, and 
harshness (NVH) during operation. Further, the invention accurately models the torque 
characteristics of the electromagnetic retarders, including through an improved saturation 
torque estimation technique, that again improves the control and operation of the overall 
system. 

[0073] The foregoing discussion discloses and describes an exemplary embodiment of the 
present invention. One skilled in the art will readily recognize from such discussion, and from 
the accompanying drawings and claims that various changes, modifications and variations can 
be made therein without departing from the true spirit and fair scope of the invention as defined 
by the following claims. 
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